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Abstract. We present A'-band commissioning observations of the Mira star prototype o Cet obtained at the ESO Very Large 
Telescope Interferometer (VUTI) with the VINCI instrument and two siderostats. The observations were carried out between 
2001 October and December, in 2002 January and December, and in 2003 January. Rosseland angular radii are derived from 
the measured visi bilities by fitting theoretical visibility functions obtained fr om center-to-limb intensity variations (CLVs) of 
Mira star models jBessell et al]|l99d:lHofmann et alj|l998tlTei et alJl2003bft . Using the derived Rosseland angular radii and 
the SEDs reconstructed from available photometric and spectrophotometric data, we find effective temperatures ranging from 
Teff = 3 192±200 K at phase O = 0. 13 to 2918± 183 K at O = 0.26. Comparison of these Rosseland radii, effective temperatures, 
and the shape of the observed visibility functions with model predictions suggests that o Cet is a fundamental mode pulsator. 
Furthermore, we investigated the variation of visibility function and diameter with phase. The Rosseland angular diameter of 
o Cet increased from 28.9 ± 0.3 mas (corresponding to a Rosseland radius of 332 ± 38 Rq for a distance of D = 107 ± 12 pc) at 
<1> = 0.13 to 34.9 ± 0.4 mas (402 ± 46 Rq) at O = 0.4. The error of the Rosseland linear radius almost entirely results from the 
error of the parallax, since the error of the angular diameter is only approximately 1%. 

Key words, instrumentation: interferometers - techniques: interferometric - stars: late-type - stars: AGE and post-AGB - 
stars: fundamental parameters - stars: individual: Mira 



1. Introduction 

Mira stars are long-period variables (LPVs) which evolve along 
the asymptotic giant branch (AGB), and are characterized by 
stellar pulsation with amplitudes as large as AV ~ 9 and well- 
defined pulsation periods (80-1000 days). In recent years, the 
comparison of theoretical pulsation models with MACHO ob- 
servations of LPVs in the LMC, in particular the reproduction 
of period ratios in multimode pulsators, has shown that Miras 
are fundamental-mode pulsators ( Wood et al, 1999) . However, 
radius measurements of Mira variables when compared to the- 
oretical pulsation calculations have generally yielded the large 
values expected for first overtone pulsators (e.g.. iFeastiil996; 
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Ivan Belle et aljfl996l) . There is clearly a problem with the in- 
terpretation of radius measurements that needs examination. 
High-resolution interferometric studies of Mira stars allow the 
determination of the size of the stellar disk, its center-to-limb 
intensity variation, surface inhomogeneities, and the depen- 
denc e of diamete r on wavelength and variabi l ity phase (see 
e.g., 'Pease' 'l93l'; 'Bonneau & Labevrie' 'l973'; Labevrie et al 
1^977; Bonneau etal. 1982; Karovska et al. 1991; Haniff etal 
1992 
11994 
U996 
1999 

i2000t iThompsonet alJl2002t .MeisnetlbOOl) . The results of 
such interferometric observations can be compared with pre- 
dictions from theoretical models of stellar pulsation and the 
atmosphere of Mira stars (e.g.. 'Watanabe & Kodairal Il979l; 
pcholz 1985; Bessell et al. 1989; Bessell et al. 1996 = BSW; 
Hof mann et al] (l99§ = HSW; fl^i et ali l2003bl = TLSW; 
llreland et al.ll2004l = LSW). Confrontation of detailed theoret- 
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ical models with high-resolution observations is crucial for 



improving our understanding of the physical properties of 


Mira stars (e.g., Hofmann et al. 


2000a; Hofmann et al. 


2001; 


IWeigeh et alJ2003: Scholz 2003 


). 



In this paper we present ESO VLTI/VINCI visibility measure- 
ments of o Cet and compare the measured visibility shape and 
the phase dependence of the visibility with model predictions. 
o Cet, the prototype of oxygen-rich Mira stars, is a very suit- 
able object for these studies, since VINCI observations exist for 
different phases and baselines, its distanc e is known (revise d 
HIPPARCOS distance 107.06 + 12.26 pc.'KnaDD e t al]l2003h . 
and a large amount of spectroscopic and photometric data is 
available for different phases. 

2. Observations and data reduction 

2.1. Observations 

A total of 48 visibility measurements of o Cet were carried out 
with VINCI (Kervella et al. 2003) at the VLTI in the commis- 
sioning period between 2001 October and 2003 January. All 
this data has been publicly released. Projected baselines rang- 
ing from 7.5 to 16.0 m were employed. VINCI is a fiber-optics 
beam combi ner instrument based on the concept of the FLUOR 
instru ment (iCoude du Foresto & Ridgwavll992tlMariotti et alJ 
With the single-mode fibers to spatially filter the wave- 
fronts perturbed by atmospheric turbulence, this beam com- 
biner provides accurate visibility measurements in spite of 
time-variable atmosphere conditions. 

Table [2 gives an overview of the VINCI observations. 
Figure[3shows the visual light curve of o Cet together with the 
dates at which the VINCI observations were caiTied out ( Matt^ 
l2003h . Figure 121 shows the visibilities of o Cet vs. spatial fre- 
quency measured at six different phases. The comparison be- 
tween the observations and the different Mira star model series 
will be discussed in Sect.|3l 

In addition to the VINCI observations, we recorded speckle 
interferograms of o Cet on 2003 October 7 with the SAO 6 m 
telescope in Russia. The speckle camera used for the observa- 
tions was equipped with a Rockwell HAWAII array. The field 
of view of the recorded speckle interferograms was 1 r.'02. A 
filter with a central wavelength of 2086 nm and a width of 
20 nm was used. The exposure time per frame was 10 ms. The 
complete data set consists of 100 speckle interferograms of o 
Cet and 360 of an unresolved reference star (HD 14652). The 
seeing (FWHM) was ~ 3'.'4. The visibilities w ere obta ined us- 
ing the speckle interferometry method ( Labevriell970l) . 

2.2. Data reduction 

Processing of the raw VINCI data for estimating the coherence 
factors jj (i.e. the fringe contrast or uncalibrated visibility) was 
caiTied out with the VINCI data reduction software provided 
by the European Southern Observatory' (version 2.0.6), based 
on wavelets transforms to derive the power spectral density. 
The processing of a series of OPD scans yields the uncali- 
brated squared visibility yu^ for each individual scan, and the 

' www.eso.org/~pballest/vinci 



Table 1. Summary of VINCI commissioning observations of o 
Cet: Date, Julian Date JD, cycle and visual phase O (see Fig.0, 
number of visibility data points A^, projected baseline length 
Bp, and baseline projection angle P.A. 



Date 


JD 


cycle+O 


N 


Sp 


[m] 


P.A. ["] 


2001 Oct 22 


2452205 


0.13 


13 


11.5 


- 16.0 


62-73 


2001 Oct 23 


2452206 


0.13 


12 


11.5 


- 16.0 


63-73 


2001 Nov 09 


2452223 


0.18 


6 


14.5 


- 16.0 


71-73 


2001 Nov 15 


2452229 


0.18 


3 


15.5 


- 16.0 


72-73 


2001 Nov 17 


2452231 


0.18 


3 


14.5 


- 15.0 


70 


2001 Dec 06 


2452250 


0.26 


6 


15.0 


- 16.0 


72-73 


2002 Jan 20 


2452295 


0.40 


2 


14.0 


- 14.5 


72 


2002 Dec 20 


2452629 


1.40 


2 


7.5 


-8.0 


73-74 


2003 Jan 09 


2452649 


1.47 


1 


g 


.0 


73 
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Julian Date -2452100 

Fig. 1. Visual light curve for o Cet llMatteill2003l) between JD 
2452100 and 2452700. The VINCI visibiHties were measured 
at the following phases: O = 0.13 (2001 Oct 22/23), O = 0.18 
(2001 Nov 09/15/17), <D = 0.26 (2001 Dec 06), 0) = 0.40 (2002 
Jan 20), O = 1.40 (2002 Dec 20), and <D = 1.47 (2003 Jan 09). 
Phases larger than 1 indicate that the observations belong to the 
next pulsation cycle. 

averaged value (i/) is used to derive the calibrated visibility 
dCoude du Foresto et al.ll997tlPerrinl2'003l) . 

The squared transfer functions T^, which account for instru- 
mental and atmospheric effects, were derived from measure- 
ments of calibrator stars and their known angular diameters. 
The expected calibrator visibility y^ai was calculated from uni- 
form disk (UD) angul ar diameters of the calibrators given in 
iRichichi & PercheronI SToOlii . and the transfer function was 
evaluated as 

y,2 ^ ^ 

y2 ' 
cal 

where jUcai is the measured uncalibrated visibility of the calibra- 
tor. The calibrators of the nights on which o Cet was measured 
are listed in Table|2 According to Eq. Q, the eiTors of the indi- 
vidual transfer function values were calculated from the errors 
of the raw visibilities and the errors of the uniform-disk visi- 
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Fig. 2. Visibilities of o Cet vs. spatial frequency measured at phase 0.13 (top left, observation dates 2001 Oct 22/23), at phase 
0.18 (top right, observation dates 2001 Nov 09/15/17), at phase 0.26 (middle left, observation date 2001 Dec 06) together 
with SAO speckle interferometry data measured at the same phase but at a different cycle (observation date 2003 Oct 07; data 
not used for the model visibility fit), at phase 0.40 (middle right, observation date 2002 Jan 20), at phase 1 .40 (bottom left, 
observation date 2002 Dec 20), and at phase 1.47 (bottom right, observation date 2003 Jan 09). The insets show an enlargement 
of the relevant spatial frequency range. Here, the error bars are included. The fits with different Mira star models are discussed in 
Sect.|3 Solid lines represent the P models from TLSW and ISW, dashed lines the D models from BSW, dotted lines the M models 
from HSW, and the dash-dotted lines represent the simple uniform-disk model CLVs. Note that from our SAO observations there 
is no indication for a visibility contribution from a circumstellar dust shell at low spatial frequencies. 
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Fig. 3. Linearly fitted squared transfer functions (see text) as a function of time (MJD = Modified Julian Date) obtained from 
calibrator stars measured (see Table|2}. For comparison, in the nights of 2001 Dec 06, 2002 Jan 20 and 2003 Jan 09 we included 
a constant fit of the squared transfer functions. The arrows mark the observation times of o Cet. 



Table 3. Properties of Mira model series from BSW, HSW, ISW, and TLSW (see text). L is the luminosity, M the stellar mass, 
and P the pulsation period in days. Rp is the Rosseland radius of the initial non-pulsating "parent" model of the Mira variable, and 
a is th e mixing-length parameter of convection (=the mixing-length in units of the local pressure scale height, see .Bohm-Vitensei 
Il958l). 



Series 


Mode 


P(d) 






Rp/Ro 


Teff/K 


a 


D 


f 


330 


I.O 


3470 


236 


2900 


1.76 


E 





328 


1.0 


6310 


366 


2700 


1.26 


P 


f 


332 


1.0 


3470 


241 


2860 


2.06 


M 


f 


332 


1.2 


3470 


260 


2750 


1.73 








320 


2.0 


5830 


503 


2250 


0.93 



bilities of the calibrators which arise from the uncertainties of 
the corresponding uniform-disk diameters. For most of the in- 
dividual transfer function values, the eiTor contributions of the 
raw visibility and the uniform-disk visibility of the calibrator 
are of the same order, resulting in total eiTors of typically ^ 5 
percent (see Fig.l^. Only in a few cases the error contribution 
of either the raw visibility or the uniform-disk visibility of the 
calibrator is significantly larger, resulting in much larger errors 



of the transfer function values. 

In order to obtain the squared transfer function for the time 
of the observation of the object, the measured values of the 
squared transfer function of each night were plotted against 
the time of observation and fitted by a straight line. Since the 
transfer function proved to be a slowly varying function which 
was rapidly sampled, as Fig.l^shows, we did not apply higher 
order polynomial fits. The error of the time-dependent transfer 
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function is given by the error of the Unear interpolation. To es- 
timate the error introduced by our method of determining the 
transfer function at the time of the observation of o Cet, we also 
applied constant fits to the three nights where the variation of 
the transfer values was largest (2001 Dec 06, 2002 Jan 20, and 
2003 Jan 09). The difference obtained for the o Cet calibrated 
visibilities using the linear and the constant transfer function fit 
procedures is 1.5% at most, leading to a final error of 0.1 mas 
for the Rosseland diameter of o Ceti. 

Since VINCI uses a K broad-band filter, model visibilities cal- 
culated with VINCI's filter response function would be ap- 
propriate for the comparison of models and measurements. 
However, those polychromatic VINCI-filter models cannot be 
calculated since monochromatic Mira star CLV models are 
not readily available. Therefore, instead of using VINCI-filter 
model CLVs, we adopted the following approximative scheme: 
In our visibility derivations we assumed the same effective 
wavelength A^s - 2.2 jjm for all calibrators and science ob- 
jects. This effective wavelength results from the assumption 
of a simple rectangular response function for the /T-band filter 
used (central wavelength 2.2 //m and bandwidth 0.4 yum) and a 
constant spectrum of the calibrators/science objects within the 
K band. This simplification introduces only a minor error in 
our calculations: the effective wavelength A^s will be different 
from 2.2 yum since t he filter response f unction is not rectangular 
in shape (see, e.g.. lWittkowski etal 1 12004 ). and the calibrator 
stars have diverse spectra. Our calculations using the calibrator 
star with the earliest spectral type (Sirius) lead to an error of 
the calibrated object visibility of 0.23% at 16 m baseline (the 
largest baseline for the observations discussed in this paper). 
The final errors (Fig. |2j of the calibrated visibilities were de- 
rived from simple Gaussian error propagation of the error of 
the fits to the measured transfer function values and the error 
of the uncalibrated visibility together with the additional effec- 
tive wavelength error 

3. Comparison of the observations witli IVIira star 
models 

3.1. Mira star models 

All Mira star models used in this paper are from BSW (D and 
E series), HSW (P, M, and O series), TLSW (P series), and 
from ISW (P series). They were developed as possible repre- 
sentations of o Cet and hence have periods P close to 332 days. 
These models dififer in pulsation mode, stellar mass M, parent 
star radius Rp (i.e., radius of the initial non-pulsating model), 
and luminosity L. The parent-star luminosities of the model se- 
ries had been chosen based on considerations of the period- 
luminosity-relation for the LMC and a theoretical metallicity 
correction (BSW, HSW). Solar abundances were assumed for 
all models. The five model series represent stars pulsating in 
the fundamental mode (f; D, P, and M) or in the first over- 
tone mode (o; E and O). Table|3]lists the properties of the non- 
pulsating parent stars of these Mira models, and Table0]shows 
the Rosseland radius and effective temperature associated with 
all P-series models. All models are named after the cycle (see 
Fig-fflfor the definition of cycle) and phase they represent. For 



Table 2. List of calibrators used (see iRichichi & PercheronI 
i2002lD and their UD angular diameters d^^. 



calibrators 


rfj^P [mas] 


nights 


sp. type 


30 Psc 


7.20 ± 0.70 


2001 Oct 22 


M3III 


T Aqr 


4.97 ± 0.05 


2001 Oct 22,23 


K5III 


y Eri 


8.51 ±0.09 


2001 Oct 22,23, 


Mllllb 






2001 Nov 09,15, 








2001 Dec 06, 








2002 Dec 20 




a Cet 


11.52 ±0.29 


2001 Oct 22, 


M1.5IIIa 






2002 Dec 20 






6.70 ±0.15 


2001 Oct 23 


M3III 


T Pup 


4.49 ± 0.07 


2001 Oct 23, 


Kim 






2002 Jan 20, 








2002 Dec 20, 








2003 Jan 09 




n Eri 


4.80 ±0.50 


2001 Nov 09,15,17 


Ml III 


Sirius 


5.60 ±0.15 


2001 Nov 09,15,17, 


AlV 






2002 Jan 20, 








2002 Dec 20 




ySCet 


5.18 ±0.06 


2001 Oct 22,23, 


KOIII 






2001 Nov 15,17, 








2001 Dec 06, 








2002 Dec 20 




<5CMa 


3.29 ± 0.46 


2001 Nov 15 


F8Iab 


HD 12274 


5.30 ±0.50 


2002 Jan 20 


K5/M0III 


eCMa 


4.13 ±0.40 


2002 Jan 20, 


K4III 






2002 Dec 20 




a Hya 


9.44 ± 0.90 


2002 Jan 20, 


K3II-III 






2003 Jan 09 




n Leo 


4.78 ± 0.26 


2002 Jan 20 


M2III 


r] Cet 


3.35 ± 0.04 


2002 Dec 20 


K1.5III 


£ Lep 


5.90 ± 0.06 


2001 Oct 22,23 


K4III 






2002 Dec 20, 








2003 Jan 09 




f Hya 


3.10 ±0.20 


2003 Jan 09 


G9II-III 


( Hya 


3.41 ± 0.05 


2003 Jan 09 


K2.5III 


yu Hya 


4.69 ± 0.50 


2003 Jan 09 


K4III 


ySCnc 


4.88 ± 0.03 


2003 Jan 09 


K4III 



example, P28 denotes the P model calculated near visual phase 
0.8 in the second cycle. The suffix "n", (e.g., PI In), denotes the 
new models from the recent fine phase grid of ISW. 
In this paper we use the conventional photospheric stellar 
Rosseland radius Rross given by the distance from the star's 
centre at which the Rosseland optical depth equals unity. This 
Rosseland radi us is no observable quan tity (cf. the discus- 
sion in Rasche k etalJll99lt ISchol7ll2003l) . Usually, it is close 
to the near-infrared radius measured in near-continuum band- 
passes except at near- minimum phases with strong molecular 
bandpass contamination (HSW, Ijacob & Schol jl2002l ISchold 
I2OO3I) . Therefore, for better comparison Table|4]also gives the 
RiM radius that is defined by the position of the t-1 layer in 
the 1 .04 fim continuum window. 

In order to derive the model visibilites, we averaged the TiT-band 
model CLV over all monochromatic CLVs between 2.0fj.m 
and 2.4 yum. An optimal comparison with VINCI visibilities, 
however, would require a slightly different calculation of the 
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Table 4. The variability phase Rosseland radius Rross in 
units of the parent star radius Rp, A'-band radius Rk in units of 
the Rosseland radius, radius /?i.o4, which is defined by the posi- 
tion of the T= 1 layer in the 1 .04 fj.m continuum window in units 
of Rp, the effective temperature Tf.ff{R), and the luminosity L for 
all P-series models from HSW and TLSW. Models denoted by 
"n" are from ISW with finer phase spacing (see text). Phases in 
the 1+0.0 to 3+0.0 range were re-assigned as described in ISW. 



Model 


cycle+O 


"Ross/^p 


n in 
"/r/"Ross 






L 












[K] 




P05 


0+0.5 


1.20 


0.65 


0.90 


2160 


1650 


P08 


0+0.8 


0.74 


0.97 


0.74 


3500 


4260 


PIO 


I +0.00 


1.03 


0.94 


1.04 


3130 


5300 


PI In 


I+O.IO 


1.17 


0.88 


1. 19 


2990 


5650 


PI2 


1+0.23 


1.38 


0.78 


1.30 


2610 


4540 


PI3n 


1+0.30 


1.53 


0.69 


1.26 


2310 


3450 


PI An 


1 -1-0 AO 


1 7^ 


U.JO 


1 1 Q 




zyzu 


PI5n 


1+0.50 


1.88 


0.49 


0.84 


1800 


1600 


PIS 


1+0.60 


1.49 


0.47 


0.85 


1930 


1910 


PIS 


1+0.87 


0.77 


0.94 


0.77 


3520 


4770 


P20 


1+0.99 


1.04 


0.92 


1.04 


3060 


4960 


P21n 


2+0.11 


1.23 


0.87 


1.21 


2790 


4750 


P22 


2+0.18 


1.32 


0.82 


1.26 


2640 


4400 


P23n 


2+0.30 


1.36 


0.78 


1.24 


2470 


3570 


P24n 


2+0.40 


1.38 


0.72 


1. 16 


2210 


2380 


P25 


2+0.53 


1.17 


0.66 


0.91 


2200 


1680 


P28 


2+0.83 


0.79 


0.97 


0.79 


3550 


5200 


P30 


2+0.98 


1.13 


0.89 


1. 14 


3060 


5840 


P35 


3+0.5 


1.13 


0.63 


0.81 


2270 


1760 


P38 


3+0.8 


0.78 


0.96 


0.78 


3570 


5110 


P40 


4+0.0 


1.17 


0.90 


1. 16 


2870 


4820 



3.2. Comparison of observed and model Rosseland 
linear radii 

Linear Rosseland radii of o Cet were obtained using the de- 
rived angular radii together with Mira's revised HIPPARC OS 
parallax of 9.34 + 1.07 mas dKnann et al.l2003UF.SAll 997l) . 
Figure|4]shows the comparison between the derived Rosseland 
linear radii and those of the corresponding models which 
provided the CLVs. The x axis represents different combina- 
tions of model series, phase, and cycles. Clearly, the E and 
O first-overtone series and the D fundamental-mode series 
predict radii that are systematically either too small (D series) 
or too large (O and E series). The fundamental-mode M series 
shows good agreement for the phases O = 0.13 - 0.26 and 
1.47, whereas the predicted radii are systematically too small 
at the phases 0.4 and 1.4. The disagreement can be partially 
explained by the lack of M models for phases close to 0.4, 
since we use near-minimum (M05 and Ml 5) models for the 
comparison. The absence of a fine phase grid also applies to 
the D and O model series, but agreement would probably not 
be improved with a finer phase grid, considering the large 
discrepancy between near-phase radii predictions and derived 
radii for O =0.18, 0.26, and 1.47. The P-series models from 
HSW, TLSW, and ISW predict, in most cases (i.e., 10 out of 
14 cycle-phase combinations), Rosseland linear radii which 
agree with the measured values within the error bars. This 
good agreement of the P models is independent of the phase 
grid spacing. The Rosseland linear radii predicted from the 
M models agree with the measured values (within the error 
bars) in 8 out of 12 cycle-phase combinations. This radius 
comparison favors the fundamental-mode pulsation P and M 
models as possible representations of Mira, whereas all first 
overtone pulsation model series can be ruled out as models 
that reliably predict the Rosseland radius of o Cet. 



polychromatic model visibility function: for each monochro- 
matic model CLV within the K band, first the squared visibility 
function has to be calculated, and then averaging over the en- 
tire wavelength range will deliver the correct squared visibility 
function. The square root of this average corresponds to t he vis- 
ibiUty curve measured with VINCI (see, e.g..lwittkowski ~t alJ 
12004 *). Comparison of the two methods for obtaining model 
visibility curves shows only small discrepancies in the relevant 
spatial frequency range and yields a systematic visibility differ- 
ence of less than 0.5%, which we included in the error budget 
of the calibrated visibility. 

The Rosseland angular radii R^^^^ of o Cet were determined 
in the following way: first, we derived the /T-band angular ra- 
dius 7?^, which is defined as the FWHM of the model CLV. 
The R^f, values (corresponding to different model-phase combi- 
nations) were determined by least-square fits between the mea- 
sured Mira visibilities and the visibilities of the correspond- 
ing model CLVs. Then, the Rosseland angular radii R^^^^ were 
derived from the obtained /iT-band angular radii and from the 
theoretical ratios RrosJRk provided by the models. For a more 
detailed description of the visibility fitting procedure we refer 
to HSW. 



3.3. Determination of T^ff 

The effective temperature Teft of a star can be derived from 
its Rosseland angular diameter and its bolometric flux (cf. 
Scholz & Takeda 1987; Bascheketal. 1991; Scholz 2003 for 
the issue of T^ff definition). Using convenient units, the rela- 
tion between T^ff, Fboi, and Rosseland angular diameter d^^^^ 
can be written as follows: 

r,ff = 2341 ( ^^^^1' ( ^] ' K. (2) 

\lO-iiWm-2/ \mas/ ^ ^ 

In order to estimate the bolometric flux, wejised photomet- 
ric measurements of o Cet in the visual (Mattel 2003) and, 
when available, in the NIR in combination with the ISO spec- 
trum ranging from 2.36 fim to 45.38 //m. In the region be- 
tween 0.78 12 yum and 1 .0834 ^m spectrophotometric data were 
taken from measurements of Wing. ( ,2003J. The 1 .25 to 2.2 jum 
data were taken from JHK photometry bv lYudinlll2003l) . Since 
the ISO spectrum was taken at maximum light ((t = 0), we 
used light curves (.Smith et alJ>200 2) obtained wi th the Diffuse 
Infrared Background Experiment (DIRBE, see ISmithl l2003l) 
to estimate the variability of o Cet in the ISO spectral range. 
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Fig. 4. Derived Rosseland linear radii for all epochs of observation and the Rosseland linear radii predicted by the nearest-phase 
models from all series. The Rosseland linear radii are given for the six measured phases O = 0.13 (top left), O - 0.18 (top 
right), (b = 0.26 (middle left), O - 0.40 (middle right), O = 1.40 (bottom left), and = 1.47 (bottom right). Open squares 
show the model values, while the filled squares with the error bars give the values derived from the observations. The large error 
bars of 10 - 15% are mainly due to the parallax error because the error contribution of the interfero metric measurements is only 
approximately 1%. 



DIRBE light curves are available in seven infrared bands be- 
tween 1.25 fim and 60 fim, giving us important complementary 
photometric data. Additional data used and the overall shape of 
the spectral energy distributions (SEDs) can be seen in Fig.|5l 
The light curves in the ISO range show that flux variations 
in the mid- and far-infrared between maximum and minimum 
light is small but not negligible. Since most of the bolomet- 
ric flux comes from the NIR, the error of Fboi induced by flux 
density changes in the ISO range makes a contribution of only 
~ 5 - 10% to the total error. Because of the small interstellar ex- 
tinction in the V and bands (Ay = 0.09, .Robertson & Feast 
Il98ll and = 0.05. iKnapp et alJl2003h and the proximity 
of Mira, no correction for interstellar extinction was applied. 
Integration of the SED was made by simple trapezoid inte- 
gration and yields Fboi = 2.63 ■ 10 Wm ^ at <D = 0.13, 
Fboi = 2.41 ■ 10-** Wm-2 at (D = 0.18, and Ffi = 2.36 ■ 10"^ 
Wm~^ at O = 0.26. The total error of the bolometric flux was 
estimated to be 25%. The error sources are mainly the cycle 
and phase uncertainties of the spectrophotometric and photo- 
metric observations. Since spectrophotometric and photomet- 
ric data were not sufficiently available for reconstructing SEDs 
at phases later than 0.3, we derived T^ff from the Rosseland an- 
gular diameter d^^_^^ and the SED according to Eq. (|2j only for 
the phases O - 0.13, 0.18, and 0.26. The error in T^ff is given 
by 



AT, 



eff 



! eff 



Adf, 



•'Ross 



1 



AF, 



bol 



(3) 



and is 6% when errors of approximately 25% and 1% are as- 
sumed for Fho\ and d^^^^, respectively. 

Figure |6l presents T^ff values derived at phases close to 
the VINCI observations from the bolometric flux and the 
Rosseland angular diameters. The T^ff values of the D model 
series are very close (within the error bars) to those derived 
from the observations at all phases. Although the PI In model 
shows agreement between observation and theory for phase 
0.13, all other cycle-phase combination of the P-model series 
as well as all other model series yield large differences between 
theoretical and measured effective temperature values. 
We also determined the luminosity range of o Cet for the 
observed phases from its bolometric flux and HIPPARCOS 
parallax (9.34 + 1.07 mas. Knapp et al. 2003) . We obtained 
936O + 314OL0forO = 0.13 and 8400 + 2820 Lq for O = 0.26. 
These luminosities are higher than those predicted by all mod- 
els for the respective phases except for the E model series, for 
which a higher luminosity of the parent star was adopted (see 
TableE). 

3.4. CLVs and visibility shapes 



As discussed in Sect. 13.21 there is good agreement between the 
derived Rosseland linear radii and the Rosseland linear radii 
predicted by model series P and M series models, while there is 
little or no agreement for the other models. On the other hand, 
the effective temperatures derived from observations are close 
to those predicted by the D-model series. In this section we in- 
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Fig. 5. Spectral energy distributions (SEDs) of o Cet for phase 
(D = 0.13 (top), O = 0.18 (middle), and O = 0.26 ^bottom ). 
The panels show visual flux densities from Mattel (2003) at 
the corresponding phase and cycle (triangles), together with the 
ISO spectrum (longward of 2.36 fj.m, dotted Une) taken at max- 
imum light. Also included are JHK photometry measurements 
from . Yudin ( 2003) (filled diamonds) taken at the phase and cy- 
cle which correspond to the epoch of the VINCI observ ations. 
The top and middle panels show spectrophotometry from Wins 



(|2003s) at 26 wavelengths between 0.78 jum and 1.10 //m (solid 
line) for the respective phases, and the lower panel shows pho- 
tometry data at 7 wavelengths: data at 0.44 //m and 0.70 yum are 
from Celis S, ( , 1982) ( filled circles), and data betw een 0.78 yum 
to 1.05 /im are from iLockwood & Wind ( Il97lh (soUd line). 
Further measurements at 0.70 iim and 0.90 //m for <1) = 0.23 
were taken from Barnes (1973) (open diamonds). In the top 
panel we also plotted DIRBE flux de nsities taken at max imum 
and minimum light (open diamonds) ("Smith et al."2002^ in the 
2.2, 3.5, 4.9, 12, and 25 jum photometric bandpasses. 
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Fig. 6. Comparison of the effective temperatures derived from 
the observations using Eq. (|2} at O = 0.13,0.18, and 0.26 
(filled squares) and the model effective temperatures (open 
squares). The measured effective temperatures were derived 
from the measured Rosseland angular radii (see Sect. 13. 2> and 
from our reconstructed SEDs. 



vestigate whether one of the model series can be favored by 
comparing the observed and theoretical visibility shapes. 
The measured visibilities at six different phases within the first 
and second cycle of the observations are shown in Fig. |2l to- 
gether with fitted model visibility curves for the P, M, and D 
models. The P-series fundamental pulsator mode model P21n 
shows much better agreement with the observation at O = 0.13 
than the D and M models or simple UD model CLVs. This can 
be clearly seen from the inset in Fig.|2l and a comparison of the 
reduced values for the different model fits. For model P21n 
we derived Xp2in ~ while we obtained x\)22 ~ 

41.5 and 

Xmii ~ -^^-^ ^'^^ models D22 and M22, respectively. For 
the simple uniform-disk model, which is also shown in Fig.|2l 
we found X\}Y) ~ 56.2. For all other phases shown in Fig. |2l 
the small number of visibility points and the insufficient cov- 
erage of the spatial frequency domain make it impossible to 
draw conclusions about which model series should be favored. 
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It should be noted that the agreement of the observed CLV 
shape with the P-series model CLV shape is also better than 
for the E and O first-overtone model series (reduced values 
of 22.9 and 45.4, respectively), which were akeady ruled out 
from radius and temperature comparisons. 
The visibility can also be affected by a circumstellar dust shell 
which is not included in the model calculations. A dust shell 
could cause a narrow visibility peak at low spatial frequen- 
cies. Therefore, we measured the visibilities of o Cet using our 
speckle camera at the SAO 6 m telescope. These additional vis- 
ibilities cover the low spatial frequency domain, showing that 
the use of dust-free CLV models is justified (see Fig.|2l middle 
left). 

Although the measured visibility shape for <[> = 0.13 is in good 
agreement with the predictions made by the P21n model, we 
nevertheless wish to point out that the observed visibility shape 
can also be affected by deviations from circular symmetry of 
the intensity distribution. Therefore, we studied possible effects 
of such asymmetries on the visibility shapes by assuming ellip- 
tical uniform-disk intensity distributions with ratios of minor to 
major axis between 0.7 and 1 .0. It turned out that the observed 
visibilities cannot be reproduced by such elliptical uniform disk 
intensity distributions. However, we note that this study cannot 
exclude the presence of other types of asymmetries in Mira, 
such as surface inhom ogeneities produced by giant convective 
cells as suggested b v iSchwarzschildl d 1 975l) . 



3.5. Diameter-phase relationship 

Taking the Rosseland diameters derived with the P model CLVs 
at the four different phases within the first cycle, we obtained 
the diameter-phase relationship presented in Fig.0 We do not 
consider the data within the second cycle of the observations, 
since the visibilities at <1) = L40 and <i> - 1.47 lie in a com- 
pletely different spatial frequency range (or baseline range, re- 
spectively) compared to the data for O < 1, and the derived 
diameters are very sensitive to small differences between the 
assumed model CLV shape and the actual CLV shape of the 
object. Figure clearly shows a monotonic increase of the 
Rosseland angular diameter with decreasing brightness from 
phase O - 0.13 {d'^^^^ - 28.9 + 0.3 mas, corresponding to 
a Rosseland linear radius of Rross - 332 + 38 Rq) to phase 
(D = 0.40 (t/^^^^ = 34.9 + 0.4 mas, Rt^^s. = 402 + 46 Ro). 
Therefore, as the visual intensity decreases, the /T-band diam- 
eter of o Cet increases. From our analysis we find an 18% 
increase of the diame ters o f o Ce t between O = 0.13 and 
15 = 0.40. Thompson et alJ ll2002l) found a similar diameter 
increase for the oxygen-rich Mira star S Lac. They obtained 
an increase of the angular size of 22% between maximum 
(O = 0.0) and minimum light (<1) - 0.5). This corresponds 
to a variation of approximately 14% between O = 0.13 and 
<i> = 0.40 (linear interpolation), which is close to our result for 
o Cet. 

Our results of the diameter-phase variation are summarized in 
Table |5] For comparison, UD angular diameters were also in- 
cluded. As can be seen, di, is similar to t/f,„, but the P models 

' Ross UD' 

agree much better with the measured visibility shapes than the 
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Fig. 7. Top: Linear Rosseland radii of o Cet derived from P21n, 
P22, P23n, and P24n series models (filled circles) at four dif- 
ferent phases (O = 0.13: Rros, = 332 + 38 Rq, (D = 0.18: 
iJRoss = 354 + 41 R^; O = 0.26: /Jross = 370 + 42 R^; 
O = 0.40: Rross - 402 + 46 Rq) and theoretical P2x models 
(filled squares). The diagram clearly shows a monotonic in- 
crease of radius with decreasing brightness, in line with the 
predictions from TLSW. Bottom: Derived (filled circles) and 
theoretical (filled squares) Rim radii of o Cet. 



simple UD models, as we discussed in the previous section. 
The intensity profiles of the best fitting models for the four dif- 
ferent phases, which were found to be the P21n, P22, P23n, 
and P24n series models, are shown in Fig.|8] The P models are 
among those with the lowest effective temperatures and exhibit 
extended faint wings near maximum light (O = 0.0) and pro- 
trusions near minimum light (O = 0.5). The extended wing in 
the intensity profile at maximum light is mainly caused by wa- 
ter molecules, which are one of the most abundant molecules in 



the atmospheres of oxygen-rich Miras (Q ei et al. 2003a). From 
the analysis of the ISO spectrum. I Yamamura et al.. C1999) con- 
cluded that H2O layers are more extended at maximum light. 
This is in agreement with the CLV near maximum light (top left 
panel in Fig.|8j. As the visual brightness decreases from max- 
imum to minimum phase, the effective temperature decreases, 
and the extended wings are less pronounced. From O = 0.18 to 
O = 0.40, water absorption increases strongly behind the shock 
front. Th is results (i) in a pro trusion of the CLV shape as seen in 
R Leo bv'Perri n et all lll999l) ("Gaussian-type" CLV) and (ii) in 
substantial molecular-band blanketing of continuum windows 
and an increase of Ri, , wh ose position is m arked by arrows in 

Ross' I r I I J 

Fig.|8](see the discussion in lTei et alJ2003ah . Note that the the- 
oretical 7?^^^,^, as well as the theoretical TJj that approximately 
describes "continuum pulsation" (Table|3, reaches their max- 
imum values at O ~ 0.4 and O ~ 0.2, respectively (Fig. 0, 
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Table 5. Diameter-phase relationship: the columns give (from left to right): phase O, models used for fitting, Rosseland angular 
diameter d^^^^ (derived by fitting P model visibility curves from TLSW and ISW to the observed visibilities), Rosseland linear 
radius Rross derived from the observations, theoretical P model Rosseland radius /^Ross.modei, and UD angular diameter 
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Fig. 8. Theoretical CLV shapes of the ISW model P21n plotted for the /T-band angular radius corresponding to o Cet at 
cE) = 0.13 (top left), of the TLSW model P22 at O = 0.18 (top right), of the ISW model P23n at O = 0.26 (bottom left), and 
of the ISW model P24n at <[) = 0.4 (bottom right) as a function of angular radius. The intensity is normalized to 1 at the disk 
center. For comparison, the intensity profiles of BSW model D22 (top right, dashed line) and the HSW model M22 (top right, 
dotted line) at O = 0.18 are also shown. The extended wing of the CLV for P models near maximum light is caused by water 
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where the difference is caused by the increasing blanketing of 
near-continuum windows towards the minimum phase. Due to 
the lack of observations close to minimum and maximum light 
phases and the size of the error bars of the radii, there is not 
enough evidence from the VINCI measurements to confirm this 



behavior of and Note in this context that the abso- 

Ross 1.04 

lute zero-phase point is uncertain by about 0.05 to 0. 1 due to the 
irregularities of both the observed (e.g.. Whitelock et al„„2000l) 
and the model-predicted (HSW, ISW) light curves. 
The CLVs predicted by the M models have comparable protru- 
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sions, and although the M model visibility fits to the observed 
visibilities do not show good agreement, they yield Rosseland 
linear radii which agree with the theoretical Rosseland linear 
radii. Because of higher temperatures and lower densities, the 
D models do not exhibit the protrusions seen in the P and M 
models. Therefore, the slope of the visibility is not as flat as in 
the case of P and M models. For a more detailed description of 
the different model CLVs, we refer to BSW, HSW, TLSW, and 
Iroet al.. 12003a) . 

4. Conclusions 

We presented /iT-band observations of o Cet obtained with 
the VLTl and its beam combiner instrument VINCI. From 
these VINCI observations at six different epochs we derived 
Rosseland angular radii using 5 different theoretical model se- 
ries from BSW, HSW, TLSW, and ISW. Using the derived 
Rosseland angular diameter and the SEDs reconstructed from 
various photometry and spectrophotometry data, we obtained 
effective temperatures ranging from T^s - 3192 + 200 K at 
4> = 0.13 to 2918 + 183 K for O = 0.26. 
We found that there is fair agreement between the Rosseland 
linear radii derived from the observed visibilities and those pre- 
dicted by the fundamental mode pulsation model series P and 
M, while there is no agreement for other models. The nonlinear 
pulsation models all start from a static "parent" star. The par- 
ent star for the best-fitting model series, the fundamental-mode 
P series, has a radius of ~240 Rq. It is clear from Fig. 4 that 
o Cet, and indeed the nonlinear pulsation models from the P 
series, nearly all have radii of ~300-400 Rq. Thus, the effect 
of large amplitude pulsation is to expand the surface layers of 
the star so that its apparent radius is considerably larger at most 
phases than the radius it would have if static (note that this ef- 
fect is milder with respect to a continuum radius like /?i .()4 ; cf. 
Table|3]and the discussion in ISW). This expansion of the ap- 
parent radius does not greatly affect the interior of the model 
or the pulsation period. Hence, if one compares observed radii 
of Miras in (say) a radius-period diagram with the radii and 
pulsation periods of the parent stars, then one will clearly not 
find agreement. This is the reason that the Miras were thought 
for so long to be first-overtone pulsators (see, e.g.. F east 1999 : 
iTuc hman 1999, and references therein). It is only by a detailed 
comparison of interferometric angular diameter measurements 
with models of large-amplitude, pulsating atmospheres that this 
problem has been solved. 

On the other hand, the effective temperatures derived from the 
observations are very close to the effective temperatures pre- 
dicted by the D-model series but higher than those predicted by 
the P and M models. Given the definition of effective temper- 
ature L - AncrR^T^g, this tells us that the P-series models are 
too low in luminosity, consistent with the luminosities derived 
in section l33l for o Cet (compared to the model values in Table 
13. The shape of the measured visibilities for o Cet at phase 
<£ = 0.13 are best fitted with the P model series, whereas all 
other model series and simple UD models show much poorer 
agreement with the observations. Taking all this into account, 
it is clear that a higher-luminosity, fundamental-mode model 
series is required for a more accurate modeling of o Cet. 



Furthermore, we found that the observed visibility functions 
and diameters change considerably from phase 0.13 to phase 
0.40. The Rosseland angular diameter of o Cet increases from 
'^Ross ~ ^^-^ - ^-^ (corresponding to a Rosseland lin- 
ear radius of Rross = 332 + 38 Rq) at O = 0.13 to d'^^^^ = 
34.9 + 0.4 mas (Rro,, = 402 + 46 Rq) at O = 0.4. Thus"' the 
diameter of o Cet increases by 18% between these two phases, 
which is in good agreement with the approximately 14% diam- 
eter increase derived fr om linear interpolation of the results by 
iThompson et alJ ( l2002i) for the oxygen-rich Mira S Lac. 
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